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^ , Abstract 
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■ We develop a new approach to robust adaptive beamforming in the presence of signal steering vector 
' errors. Since the signal steering vector is known imprecisely, its presumed (prior) value is used to find 

a more accurate estimate of the actual steering vector, which then is used for obtaining the optimal 
beamforming weight vector. The objective for finding such an estimate of the actual signal steering vector 
is the maximization of the beamformer output power, while the constraints are the normalization condition 
and the requirement that the estimate of the steering vector does not converge to an interference steering 
\ vector. Our objective and constraints are free of any design parameters of non-unique choice. The resulting 

C^, , optimization problem is a non-convex quadratically constrained quadratic program, which is NP hard in 

general. However, for our problem we show that an efficient solution can be found using the semi-definite 

00 ' 

relaxation technique. Moreover, the strong duality holds for the proposed problem and can also be used 
2 i for n„di„ g the optinta, solution efficient,, and at low contplexit, In soute special eases, the solution cau 

J> , be even found in closed-form. Our simulation results demonstrate the superiority of the proposed method 

•i-H . 

■ over other previously developed robust adaptive beamforming methods for several frequently encountered 
& ' types of signal steering vector errors. 
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I. Introduction 

Robust adaptive beamforming design has been an intensive research topic over several decades 
due to, on one hand, its ubiquitous applicability in wireless communications, radar, sonar, micro- 
phone array speech processing, radio astronomy, medical imaging, and so on; and on the other 
hand, because of the challenges related to the practical applications manifesting themselves in the 
robustness requirements for adaptive beamformers. The main causes of performance degradation 
in adaptive beamforming are small sample size and imprecise knowledge of the desired signal 
steering vector in the situation when the desired signal components are present in the training 
data. The traditional design approaches to adaptive beamforming [HI- flU do not provide sufficient 
robustness and are not applicable in such situations. Thus, various robust adaptive beamforming 
techniques gained a significant popularity due to their practical importance Q. The most popular 
conventional robust adaptive beamforming approaches are the diagonal loading technique J6j, 
and the eigenspace-based beamforming technique ©, 10. However, it is not clear in the former 
approach how to obtain the optimal value of the diagonal loading factor, whereas the eigenspace- 
based beamforming is known to suffer from the so-called subspace swap phenomenon at low 
signal-to-noise ratios (SNRs) ifTOl and requires exact knowledge of the signal-plus-interference 
subspace dimension. 

In the last decade, several approaches to robust adaptive beamforming based on a rigorous 
modeling of the steering vector mismatches have been developed. For the case when the mismatch of 
the signal steering vector is modeled as deterministic unknown norm bounded vector, the so-called 
worst-case-based adaptive beamforming design has been proposed in [fTT]|. |[T2l|. The relationship 
of the worst-case-based design to the diagonal loading principle with adaptive diagonal loading 
factor has been explored in ifTTTl . lfT3l as well as the generalization to the ellipsoidal steering vector 
uncertainty set has been developed in |[T4|. The uncertainty of second order statistics of the desired 
signal have been also considered in lfl"5l . If the signal steering vector mismatch is modeled as random 
unknown with known Gaussian or unknown distribution, the corresponding optimization problem 
belongs to the class of stochastic programming problems, and the corresponding probabilistically- 
constrained robust adaptive beamformer has been developed in |[T6ll . ifTTll . The relationship between 
the worst-case-based and the probabilistically-constrained robust adaptive beamformers has also 
been shown in [fTTl . lfT8l . However, some design parameters such as the norm of the signal steering 
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vector mismatch in the worst-case-based design or the acceptable beamforming outage probability in 
the probabilistically-constrained design are assumed to be known in the aforementioned techniques. 

Another approach to robust adaptive beamforming which is based on estimating actual signal 
steering vector based on the knowledge of presumed steering vector has been proposed in |[T9ll , 
ll20l . The main idea of this approach is to estimate the signal steering vector so that the maximum 
of the beamformer output power is achieved, while the convergence of the steering vector estimate 
to an interference steering vector is avoided. The latter convergence can be avoided by imposing 
different constraints. The projection constraint to the space to which the desired signal belongs is 
used in EOl and the solution based on sequential quadratic programming (SQP) is developed. The 
advantage of the method of ll20ll is that no design parameters of non-unique choice are required. 
However, the disadvantage is that the complexity of SQP is rather high, which makes the method 
less attractive for practical use. It is interesting that the robust adaptive beamforming formulation 
based on desired signal steering vector estimation has been also considered in [|2T|. but the typical 
for the worst-case-based methods [fTTIl - [TT4l norm bound constraint on the steering vector mismatch 
has been used there in order to guarantee that the estimate of the desired signal steering vector 
does not converge to an interference steering vector. As a result, the method of j[2l| is one of the 
various implementations of the adaptive diagonal loading based techniques, and it is not free of 
design parameters. 

In this paper, we develop a new robust adaptive beamforming method which is free of any design 
parameters of non-unique choice. This method is based on the signal steering vector estimation via 
beamformer output power maximization under the constraint on the norm of the steering vector 
estimate and the requirement that the estimate of the steering vector does not converge to an 
interference steering vector. To satisfy the latter requirement, we develop a new constraint which is 
a convex quadratic constraint. Then the corresponding optimization problem is a non-convex (due 
to the steering vector normalization condition) homogeneous quadratically constrained quadratic 
programming (QCQP) problem. In general, QCQP problems may not have a strong duality property, 
which leads to the situation when the solution of the corresponding semidefinite programming (SDP) 
relaxation-based problem is of rank higher than one and randomization procedures have to be used 
to find an approximate rank-one solution ll23l - ll26ll . The probability that the so-obtained rank- 
one solution coincides with the exact solution is less than one [|27l . However, in the case of our 
signal steering vector estimation problem, we show that a rank-one solution can be found efficiently 
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using the SDP relaxation technique. Moreover, the strong duality holds for the proposed problem, 
which means that a rank-one solution can also be found based on the solution of the convex 
dual problem. Some special cases and interesting relationships are also considered. Our simulation 
results demonstrate the superiority of the proposed method over other previously developed robust 
adaptive beamforming techniques. 

This paper is organized as follows. Data model, beamforming formulation, and necessary back- 
ground are given in Section II. In Section III, we formulate the problem of interest. A complete 
analysis of the problem and its rank-one solution is given in Section IV. Section V overviews some 
special cases and draws interesting relations to the existing methods. Simulation results comparing 
the performance of the proposed method to the existing methods are shown in Section VI. Finally, 
Section VII presents our conclusions. 

II. System Model and Background 

Consider a linear antenna array with M omni-directional antenna elements. The narrowband 
signal received by the antenna array at the time instant k can be written as 



where s(k), i(/c), and n(k) denote the M x 1 vectors of the desired signal, interference, and 
noise, respectively. The desired signal, interference, and noise components of the received signal 
CQ> are assumed to be statistically independent to each other. The desired signal can be written as 
s(k) = s(k)a, where s(k) is the desired signal waveform and a is the steering vector associated 
with the desired signal. 

The beamformer output at the time instant k can be written as 



where w is the M x 1 complex weight (beamforming) vector of the antenna array and (-) H stands 
for the Hermitian transpose. 

Assuming that the steering vector a is known precisely, the optimal weight vector can be obtained 
by maximizing the beamformer output signal-to-noise-plus-interference ratio (SINR) 01 



x(fc) = s(k) + i(k) + n(fc) 



(1) 



y(k) 



w tf x(Jc) 



(2) 



SINR 




(3) 
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where of is the desired signal power, R i+n = E{(i(k) + n(k))(i(k) + n(k)) H } is the M x M 
interference-plus-noise covariance matrix, and stands for the statistical expectation. Since 

Rj+„ is unknown in practice, it is substituted in © by the data sample covariance matrix 

1=1 

where K is the number of training data samples which also include the desired signal component. 
The sample version of the problem of maximizing © is known as the minimum variance (MV) 
sample matrix inversion (SMI) beamforming and it is based on its conversion to the mathematically 
equivalent problem of minimizing the denominator of © under fixed numerator, that is, 

minw^Rw subject to w^a = 1. (5) 

w 

The problem © is convex and its solution can be easily found as wmv-smi = «R _1 a where 
a = l/a H R" 1 a ffl. 

In practice, the steering vector a is not known precisely and only its inaccurate estimate p, called 
hereafter as presumed steering vector, is available. Several rigorous approaches, which address the 
problem of imprecise knowledge of the desired signal steering vector have been developed in the 
last decade. First of them assumes that the actual steering vector a can be explicitly modeled as a 
sum of the presumed steering vector and a deterministic mismatch vector S, i.e., a = p + <5 lITTI . 
lfT2l . Here d is unknown but it is known that \\S\\ < e for some bound value e, where || ■ || is 
the Euclidian norm of a vector. This approach has been also generalized for ellipsoidal uncertainty 
case in addition to the aforementioned spherical uncertainty |fl4l . Assuming spherical uncertainty 
set for 6, i.e., A(e) = {a = p + d | ||<5|| < e\, the worst-case-based robust adaptive beamforming 
aims at solving the following optimization problem 

minw^Rw subject to max jw^al > 1. (6) 

In turn, the problem © is equivalent to the following convex optimization problem [fTTTl 

minw Rw subject to e||w|| < w H p - 1 (7) 

w 

which can be efficiently solved using second-order cone programming (SOCP) 0T], ll28l or nu- 
merical Lagrange multiplier techniques 031, ||29ll . 
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Another approach to robust adaptive beamforming is based on the assumption that the vector 6 
is random. Then the problem © changes to 

minw Rw subject to Pr{| w^a| > 1} > Po (8) 

w 

where Pr{-} denotes probability and p is preselected probability value. In the case of Gaussian 
distributed 6 and the case when the distribution of 6 is unknown and assumed to be the worst 
possible, it has been shown that the problem © can be approximated by the following problem 

ini, ma, 

minw^Rw subject to e||Qy 2 w|| < w"p - 1 (9) 

w 

where Qs is the covariance matrix of random mismatch vector 6, and e = y — ln(l — p ) if S is 
Gaussian distributed or e = l/n/l — Po if the distribution of <5 is unknown. Both problems © and 
© have similar SOCP structure and can be solved efficiently. 

The third approach to robust adaptive beamforming aims at estimating the steering vector a 
based on the prior given by the presumed steering vector p lfT9~l . Il20l . The estimate of the steering 
vector a is found so that the beamformer output power is maximized while the convergence of the 
estimate of a to any interference steering vector is prohibited. Indeed, the solution of © can be 
written as a function of unknown 5, that is, w(<5) = aR _1 (p + S). Using the latter expression, the 
beamformer output power can be also written as a function of d as 

P(S) = . (10) 

(p + 5)^R- 1 (p + 5) 

Thus, such estimate of S or, equivalently, such estimate of a that maximizes (flOl) will be the best 
estimate of the actual steering vector a under the constraints that the norm of a equals \J~M and 
a does not converge to any of the interference steering vectors. The latter can be guaranteed by 
requiring that 

P ± (p + 5) = P ± a = (11) 

where P^ = I — UU ff , U = [ui, u 2 , . . . , Ul], Uj, I = 1, . . . , L are the L dominant eigenvectors of 
the matrix C = J e d(9)d H (9) dO, d{6) is the steering vector associated with direction 9 and having 
the structure defined by the antenna geometry, 6 is the angular sector in which the desired signal 
is located, 6 and a stand for the estimate of the steering vector mismatch and for the estimate 
of the actual steering vector, respectively and I is the identity matrix. The resulting optimization 
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problem is non-convex, but has been solved in ||20l using SQP technique. A similar approach based 
on steering vector estimation has been also recently taken in ll30ll for the case when a is partially 
known, for example, when array is partially calibrated, that significantly simplifies the problem. The 
following interesting relationship is also worth mentioning. If the constraint ([111) is replaced by the 
constraint \\8\\ < e used in the worst-case-based beamformers, the convergence to an interference 
steering vector will also be avoided, but the problem becomes equivalent to the worst-case-based 
robust adaptive beamforming (see EH)- 

Finally, it is worth mentioning that the eigenspace-based beamformer j8]|, flU is also based on 
correcting/estimating the desired signal steering vector. Taking the presumed steering vector p as a 
prior, the eigenspace-based beamformer finds and uses the projection of p onto the sample signal- 
plus-interference subspace as a corrected estimate of the steering vector. The eigendecomposition 
of © yields 

R = EAE^ + GTG H (12) 

where the M x J + 1 matrix E and M x M — J — 1 matrix G contain the signal-plus-interference 
subspace eigenvectors of R and the noise subspace eigenvectors, respectively, while the J+l x J+l 
matrix A and M — J — lxM — J — 1 matrix T contain the eigenvalues corresponding to E and G, 
respectively. Here, J is the number of interfering signals. Then the eigenspace-based beamformer 
is given by 

w cig = R _1 a = R^EE^p = EA _1 E H p (13) 

where a = EE H p is the projection of the presumed steering vector p onto the sample signal- 
plus -interference subspace and EE H is the corresponding projection matrix. As compared to the 
beamformer of [20] based on the estimation of steering vector, the eigenspace-based beamformer 
may suffer from a high probability of subspace swaps as well as incorrect estimation of the signal- 
plus -interference subspace dimension. 

III. New beamforming Problem Formulation 

The problem of maximizing the output power (flOl) is equivalent to the problem of minimizing the 
denominator of (flOl) . The obvious constraint that must be imposed on the estimate a, is that the norm 
of a must be equal to y~M, i.e., ||a|| = yM. This normalization condition, however, does not protect 
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the estimate a from possible convergence to an interference steering vector. In order to avoid such 
convergence, we assume that the desired source is located in the angular sector 6 = [6> min , 6> max ] 
which can be obtained, for example, using low resolution direction finding methods. The angular 
sector is assumed to be distinguishable from general locations of all interfering signals. In turns, 
the sector denotes the complement of the sector 0, i.e., combines all the directions which lie 
outside of 0. Let us define the M x M matrix C as C = /g d(9)d H (9) d6. Then constraint 



for a uniquely selected value A (see Example 1 below), will force the estimate a not to converge 
to any interference steering vector with the directions within the angular sector 0. To illustrate how 
the constraint (fl4)) works, let us consider the following example. 

Example 1: Consider uniform linear array (ULA) of 10 omni-directional antenna elements spaced 
half wavelength apart from each other. Let the range of the desired signal angular locations be 
= [0°, 10°]. Fig. Q] depicts the values of the quadratic term d H (9)Cd(9) for different angles. The 
rectangular bar in the figure marks the directions within the angular sector 0. It can be observed 
from this figure that the term d H (9)Cd(9) takes the smallest values within the angular sector 0, 
where the desired signal is located, and increases outside of this sector. Therefore, if A is selected 
to be equal to the maximum value of the term d H (9)Cd(9) within the angular sector of the desired 
signal 0, the constraint (fl4l) will guarantee that the estimate of the desired signal steering vector 
will not converge to any interference steering vectors. Note also that the constraint (TBI) is an 
alternative to the constraint (fTT|) used in EDI . However, the constraint (fTT|) may result in the noise 
power magnification at low SNRs (see EOl ). while the constraint (fl4l) helps to alleviate the effect 
of the noise power magnifying at low SNRs by not collecting the noise power from the continuum 
of the out-of-sector directions 0. 

Taking into account the normalization constraint and the constraint (fl4l . the problem of estimating 
the desired signal steering vector based on the knowledge of the prior p can be formulated as the 
following optimization problem 



a ff Ca < A 



o 



(14) 



mm 



(15) 



a 



subject to 




(16) 



a ff Ca < A 



o 



(17) 
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where the prior p is used only for selecting the sector 6. Due to the equality constraint (fl6l) , 
which is a non-convex one, the QCQP problem of type (fT5l)-(fl7l) is non-convex and an NP-hard in 
general. However, as we show in the following section, an exact and simple solution specifically 
for the problem (fl5T)-(fT7T) can be found using the SDP relaxation technique and the strong duality 
theory. 

IV. Steering Vector Estimation via Semi-Definite Programming Relaxation 

QCQP problems of type (fl5T)-<fT7l) can be solved using SDP relaxation technique. The first step 
is to make sure that the problem (fl31)-(fT7l) is feasible. Fortunately, it can be easily verified that 
(fl3T> - (fT7T) is feasible if and only if A /M is greater than or equal to the smallest eigenvalue of the 
matrix C. Indeed, if the smallest eigenvalue of C is larger than A /M, then the constraint (TlVT) can 
not be satisfied for any estimate a. However, A selected as suggested in Example 1 will satisfy 
the feasibility condition that will guarantee the feasibility of the problem (fTST) — (fTvT) . 

A. Semi-Definite Programming Relaxation 

If the problem (fl3T)-(fT7T) is feasible, the equalities a^R^a = Tr (R _1 aa ff ) and a ff Ca = 
Tr (Cak H ), where Tr(-) denotes the trace of a matrix, can be used to rewrite it as the following 
optimization problem 

min Tr(R _1 aa^) (18) 

a 

subject to Tr(aa^) = M (19) 
Tr(Ckk H ) < A . (20) 
Introducing the new variable A = ak H , the problem (fT8l)-(T20l) can be casted as 

min Tr(R -1 A) (21) 

A 

subject to Tr(A) = M (22) 

Tr(CA) < A (23) 

rank(A) = 1 (24) 

where rank(-) stands for the rank of a matrix and it is guaranteed by the combination of the 
constraints (|22|) and (|24|) that A is a positive semi-definite matrix, i.e., A y 0. 
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The only non-convex constraint in the problem (l2TT)-(l24l) is the rank-one constraint (1241) while 
all other constraints and the objective are convex. Using the SDP relaxation technique, the relaxed 
problem can be obtained by dropping the non-convex rank-one constraint (1241) and replacing it by 
the semi-definiteness constraint A y 0, which otherwise is not guaranteed if (1241) is not present. 
Thus, the problem (12"TI)-(12"4l) is replaced by the following relaxed convex problem 

min Tr(R -1 A) (25) 

A 

subject to Tr(A) = M (26) 

Tr(CA) < A (27) 

A y 0. (28) 

There are two features related to the use of SDR that have to be addressed. First, it is possible, 
in general, that the original problem is infeasible, however the relaxed one is feasible. Second, the 
optimal solution of the relaxed problem is, in general, an approximation of the optimal solution 
of the original problem. Thus, it is desirable in general to estimate the approximation bounds for 
the approximate solution and the probability that both approximate and exact optimal solutions 
coincide j|2~7l . In the following we will address these issues. 

B. Feasibility and Rank of the Optimal Solution 

The result that connects the feasibility of the relaxed problem (T25T)-(l28l) to the feasibility of the 
original problem (fl"5l)-([T7T) is given in terms of the following theorem. 

Theorem 1: The problem (l25l)-(l28T) is feasible if and only if the problem (fT5l)-(fT7l) is feasible. 
Proof: See Appendix. 

If the relaxed problem (T25T)-(|28T) has a rank-one solution, then the principle eigenvector of the 
solution of (I25l)-(l28l) will be the exact solution to the original problem (fT5l)-(fT71). Otherwise, 
randomization procedures 11231 . Il2~6ll have to be used, which can find the exact optimal solution of 
the original problem only with a certain probability [|27l . However, under the condition that the 
original optimization problem (fT5l)-(fT7l) or, equivalently, the relaxed problem (12~5l)-(l28l) is feasible, 
the solution of the original problem can be extracted from the solution of the relaxed problem by 
means of the following constructive theorem. 

Theorem 2: Let A* be the rank r optimal minimizer of the relaxed problem (l25l) - (12~8l) . i.e., 
A* = YY H (where Y is an N x r matrix). If r = 1, the optimal solution of the original problem 
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simply equals Y. Otherwise, it equals Yv, where v is an r x 1 vector such that ||Yv|| = yM 
and v ff Y^CYv = Tr(Y H CY). Then one possible solution for the vector v is proportional to 
the sum of the eigenvectors of the following r x r matrix 

D = ±-Y H Y - _ . (29) 

M Tr {Y H CY) 



Proof: See Appendix. 

Finally, we can prove the following result on the uniqueness of the rank-one solution of the 
relaxed problem (I25l)-(l28l). 

Theorem 3: Under the condition that the solution of the original optimization problem (fl3T) - 
([171) is unique regardless of a phase shift, where the latter means that if a and a' are both optimal 
solutions, then there exists such phase shift <p that a = a'e- 7 ^, the solution of the relaxed problem 
ft25\)-ft28\) always has rank one. 

Proof: See Appendix. 

Note that the phase shift plays no role in the desired signal steering vector estimation problem 
since the output power (flOl) as well as the output SINR do not change if a undergo any phase 
rotation. Thus, the uniqueness condition regardless a phase shift in Theorem 3 is proper . Under 
this condition, the solution of the relaxed problem (I25l)-(l28l) is rank-one and the solution of the 
original problem (fT5l)-(fT7l) can be found as a dominant eigenvector of the optimal solution of the 
relaxed problem (T25l)-(l28l). However, even such uniqueness condition regardless a phase shift is not 
necessarily satisfied for the problem (TT3T) — (TTvT) (see Example 2 below), and then we resort to the 
constructive Theorem 2, which shows how to find the rank-one solution of (TT3T) — (TTvT) algebraically 
without any use of randomization procedures. 

Example 2: As an example of the situation when (fT5T)-(fl7T) does not have a unique solution, let 
us consider a ULA with 10 omni-directional antenna elements. The presumed direction of arrival 
of the desired user is assumed to be 9 P = 3° with no interfering sources and the range of the 
desired signal angular locations is equal to 6 = [9 P — 12°, 9 P + 12°]. The actual steering vector of 
the desired user is perturbed due to the incoherent local scattering effect and it can be expressed 
as a = p + b, where p = d(3°) is the steering vector of the direct path and b is the steering 
vector of the coherently scattered path. Let us consider the case when b is orthogonal to p. This 
later condition can be satisfied if b is selected as d(— 8.5°). In this case, both of the vectors p 
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and b are the eigenvectors of the matrix Rr 1 which correspond to the smallest eigenvalue. Since, 
these vectors satisfy the constraints (fT6l)-(fl7T) and correspond to the minimum eigenvalue, both of 
them are optimal solutions of the optimization problem (fT5l)-(fT7T), thus, the solution of the problem 
(fTSb— (fTTT) is not unique. 

C. Solution Based on Strong Duality 

The solution of (fl"5T)-([T7T) can also be found using the strong duality theory. It follows from 
Theorem 2 that the optimal value of the relaxed problem (|25T)-(|28T) is the same as the optimal value 
of the the original problem (fT5l)-(fl7l). It is because in addition to the fact that Yv is the optimal 
solution of the original problem (fT5l)-(fT7l) (see Theorem 2), Yv(Yv) fl is also the optimal solution 
of the relaxed problem (|25T)-(|28T) (see the proof of Theorem 2). Furthermore, the dual problem of 
the the relaxed problem (|25T)-(|28T) is the same as the dual problem of the original problem (fT5l)-(fT7l). 
Indeed, by maximizing the dual function of the problem (l2~5T)-(l28T). which is the same as the dual 
function of (fT5T)-([T7T). the dual problem for both the relaxed and original problems can be written 
as 

max 71M - 7 2 A (30) 

71:72 

subject to RT 1 - 7x1 + 72 C ^ (31) 

where 71 and 72 > are the Lagrange multipliers associated with constraints (fT6l) and (fTTT) of the 
original problem or the constraints d26l) and (l27l) of the relaxed problem, respectively. Since the 
relaxed problem (|25 |) - (|28|) is convex, the strong duality between (|25T> - (|28l) and <|30])-<|3D holds, 
i.e., the optimal value of (I25l)-(l28l) is the same as the optimal value of (ISljlHTJlT). It implies that 
the optimal value of the dual problem (f30b — (T3TT) is also the same as the optimal value of the 
original problem (fT5l)-([T7l). Thus, the strong duality between the dual problem (f30b— (f3Tb and the 
original problem (fT5l)-(fT7T) also holds. It is worth mentioning that the strong duality of a non-convex 
quadratic optimization problem with two positive semi-definite quadratic constraints has been also 
studied in the recent work OTTl . Particularly, it has been shown that if a non-convex quadratic 
optimization problem with two quadratic constraints is strictly feasible, then strong duality holds. 
This result agrees with our above conclusion that the strong duality between (f30l)— (f3TT) and (fl"5l)- 
(flTT) holds. Indeed, it can be easily shown that if A /M is greater than the smallest eigenvalue 
of the matrix C, then the problem (fT5l)-(fT7l) is strictly feasible, and thus the result of OTTl . Il32l 
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applies and the strong duality between (l30l)-(l3TT) and (TT3T) — dTTT) holds. Moreover, if A /M is equal 
to the smallest eigenvalue of the matrix C, the problem (fT5l)-(fl7l) has limited number of feasible 
points and it can be solved easily by checking all these points. Thus, it is simply assumed in the 
sequel that A /M is greater than the smallest eigenvalue of C. 

The dual problem (f30b— (f3Tb belongs to the class of SDP problems and, thus, can be solved 
efficiently using, for example, interior-point methods. Moreover, it contains only two optimization 
variables. Let the optimal solution of the the dual problem (I30l)-(|3TI) be 7J 1 and 7^. It is easy to 
see that 7^ is always strictly positive and the matrix Rr 1 — 7*1 + 72C is rank deficient. Indeed, in 
order to maximize the objective function (l30l) for a fixed 72 > 0, 71 should be equal to the smallest 
eigenvalue of the matrix Rr 1 + 72C which makes the matrix Rr 1 — 71 1 + 72C rank deficient. 
Furthermore, since for every nonnegative 72, Rr 1 + 72C is a positive definite matrix, we obtain 
that 7j" is positive and Rr 1 — 71I + 72C is rank deficient. Since the strong duality between (fT5l-(fT7T) 
and (l30l)-(l3TI) holds, the necessary and sufficient optimality conditions can be written as 



(R- 1 -7i*I + 7 2 *C)a = M (32) 

a^a = M (33) 

7 2 *(a^Ca- A ) = (34) 

a ff Ca < A . (35) 

where 0^/ is the vector of zeros of length M. Moreover, using the fact that 7^ is strictly positive 
and the matrix Rr 1 — 7*1 + 7|C is rank deficient, the solution of the original problem (fT5l)-(fT7l) 
can easily be foundQ There are two possible situations. 

(i) The matrix Rr 1 — 7*1 + 7^0 has only one zero eigenvalue. In this case, the only vector which 
satisfies (l32l) is given by 

a = ^/Mplir 1 -7il + 7 2 *c) (36) 



where p{-} denotes the eigenvector of a matrix which corresponds to the smallest eigenvalue. 
Indeed, (|36l) satisfies the necessary and sufficient optimality conditions (132b— (1351) and, thus, is 
the optimal solution of (IT3T) — (ITtT). 

'Note that the general form of the optimality conditions d32M35b has been solved in 1311 . 
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(ii) The matrix Rr 1 — 7JT+72C has more than one zero eigenvalue. In this case, consider the matrix 
F each column of which is an eigenvector of the matrix Rr 1 — 7*1 + 72 C corresponding to 
zero eigenvalue. Thus, the dimension of F is M x q where q is the number of zero eigenvalues 
of the matrix Rr 1 — 7*1 + 72C. Then vectors a that satisfy the condition (1321 can be written 

as 

a = Fq (37) 

where q is a k x 1 vector. The optimality conditions (I32l)-(l35l) can be rewritten then in terms 
of q and F as 

q^q = M (38) 
7 2 *(q H F ff CFq-A ) = (39) 
q^F^CFq < A . (40) 

Let fi rnax and /i min denote the largest and the smallest eigenvalues of the matrix F^CF and 
f max and tfnin stand for their corresponding eigenvectors. Then the following two subcases 
should be considered. 

(ii.a) The first subcase is when 72 = 0. Then jjL min < A /M and q can be simply chosen as 

q = y/Mf min . (41) 



(ii.b) The other subcase is when 72 > 0. Then (l39l) implies that q^F^CFq = A and as in the 
previous subcase fi rnin < A /M. Moreover, [i max > A /M. Therefore, q can be chosen 
as a linear combination of f min and i max as follows 



q = VM(VT^9f mm + V6f max ) (42) 

where 9 = (A /M - /i min )/ (fi max - l^min)- 
As soon as the estimate a is obtained, the beamforming weight vector can be straightforwardly 
computed as 

w SD p = a'R^a (43) 

where a' = l/a^R^a. The beamformer (l43l) can be compared to the eigenspace-based robust 
adaptive beamforming (fT3l) where the imprecisely known signal steering vector is corrected by 
projecting it to the signal-plus-interference subspace. However, a significant difference is that no 
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knowledge of the dimension of the signal-plus-interference subspace is needed in the proposed 
beamforming as well as no subspace swap can happen at low SNRs as in (fT3l) . 

V. Special Cases and Relationships 

A. Simple Solution Under the Constraint (fTTI) 

For high and moderate SNRs, the protection against convergence to an interference steering 
vector can be ensured by means of the constraint (fTTj) (see also Il20l0 . Then the corresponding 
desired signal steering vector estimation problem can be written as 



min & H R L k (44) 

a 

subject to ||a|| = \[~M (45) 

P x a = 0. (46) 



This problem differs from the problem (fT5l)-(fT71) only by the constraint (1461) . The problem (I44l)-(l46l) 
is a non-convex optimization problem, but a simple closed-form solution can be found. The main 
idea is to first find a set of vectors satisfying the constraint (l46l) . Note that P^a = implies that 
a = UU ff a and, therefore, we can write that 

a = Uv (47) 

where v is a L x 1 complex valued vector. Using (|47T) . the optimization problem (I44l)-(l46l) can be 
equivalently rewritten in terms of v as 

min v H U H R _1 Uv (48) 

V 

subject to || v || = \[M. (49) 

Finally, the solution to the optimization problem of type (|48l)-(|49l) is known to be the eigenvector 
of the matrix U^R -1 !! which corresponds to the minimum eigenvalue. Thus, the estimate of the 
steering vector can be obtained as 

a=U-p{u ff R _1 u} (50) 
and the corresponding beamforming vector is 

wi = a'R^U • p |u ff R _1 u} . (51) 
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As compared to the the eigenspace-based robust adaptive beamforming (fT3l) , the steering vector 
in the beamformer (|5D) is the eigenvector corresponding to the smallest eigenvalue of the matrix 
U^R^U. 

B. Signal-to-Interference Ratio —> oo 

In the case when the signal-to-interference ratio approaches infinity (SIR—)- oo), it is guaranteed 
that the estimate of the desired signal steering vector will not converge to an interference steering 
vector and, thus, the constraint (fTTT) is never active and can be dropped. Then the optimization 
problem (fl3T)-([T7T) simplifies as 

min a^R^a (52) 

a 

subject to ||a|| = \/M. (53) 

The solution of the latter problem is a = \/Mp(R _1 ). Interestingly, the optimization problem 
(I52l)-(l53l is the same as the optimization problem ET1 (39)] which is obtained after dropping the 
constraint \\S\\ < e. This relationship holds only for SIR— > oo when no additional constraints are 
required to guarantee that the estimate of the steering vector does not converge to an interference 
steering vector. In this case, the proposed and the the worst-case-based beamformer are the same, 
that is, 

w 2 = a'R _1 a = a'RrV {rt 1 } . (54) 

VI. Simulation Results 

Throughout the simulations, a ULA of 10 omni-directional antenna elements with the inter- 
element spacing of half wavelength is considered. Additive noise in antenna elements is modeled 
as spatially and temporally independent complex Gaussian noise with zero mean and unit variance. 
Two interfering sources are assumed to imping on the antenna array from the directions 30° and 
50°, while the presumed direction towards the desired signal is assumed to be 3°. In all simulation 
examples, the interference-to-noise ratio (INR) equals 30 dB and the desired signal is always present 
in the training data. For obtaining each point in the examples, 100 independent runs are used. 

The proposed SDP relaxation-based beamformer is compared with three other methods in terms 
of the output SINR. These robust adaptive beamformers are (i) the worst-case-based robust adaptive 
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beamformer ©, (ii) the SQP-based beamformer of |[20~ll , and (iii) the eigenspace-based beamformer 
(TT3T) . For the proposed beamformer and the SQP-based beamformer of EDI , the angular sector of 
interest is assumed to be 6 = [9 P — 5°, 9 P + 5°] where 9 P is the presumed direction of arrival of 
the desired signal. The CVX MATLAB toolbox [[331 is used for solving the optimization problems 
([25T)-(T28T) and (|30T>-(|3TT) and the value of A is set equal to the maximum value of the d H (9)Cd(9) 
within the angular sector of interest 6. The value 5 = 0.1 and 6 dominant eigenvectors of the matrix 
C are used in the SQP-based beamformer and the value e = 0.3M is used for the worst-case-based 
beamformer as it has been recommended in [fTTTl . The dimension of the signal-plus-interference 
subspace is assumed to be always estimated correctly for the eigenspace-based beamformer and 
equals 3. 

A. Example 1: Exactly known signal steering vector 

In the first example, we consider the case when the actual desired signal steering vector is 
known exactly. Even in this case, the presence of the signal of interest in the training data can 
substantially reduce the convergence rates of adaptive beamforming algorithms as compared to the 
signal-free training data case (8). In Fig. [2l the mean output SINRs for the aforementioned methods 
are illustrated versus the number of training snapshots for the fixed single-sensor SNR = 20 dB. 
Fig. [3] displays the mean output SINR of the same methods versus the SNR for fixed training 
data size of K = 30. It can be seen from these figures that the proposed beamforming technique 
outperforms the other techniques even in the case of exactly known signal steering vector. It is 
especially true for small sample size. 

B. Example 2: Signal Spatial Signature Mismatch Due to Wavefront Distortion. 

In the second example, we consider the situation when the signal steering vector is distorted by 
wave propagation effects in an inhomogeneous medium. Independent-increment phase distortions 
are accumulated by the components of the presumed steering vector. It is assumed that the phase 
increments remain fixed in each simulation run and are independently chosen from a Guassian 
random generator with zero mean and variance 0.04. 

The performance of the methods tested is shown versus the number of training snapshots for 
fixed single-sensor SNR= 20 dB in Fig. @]and versus the SNR for fixed training data size K = 30 
in Fig.|5l It can be seen from these figures that the proposed beamforming technique outperforms all 
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other beamforming techniques. Interestingly, it outperforms the eigenspace-based beamformer even 
at high SNR. This performance improvement compared to the eigenspace-based beamformer can be 
attributed to the fact that the knowledge of sector which includes the desired signal steering vector 
is used in the proposed beamforming technique. Fig. \5\ also illustrates the case when SNR^INR 
where INR stands for interference-to-noise ratio. This case aims to illustrate the situation when 
SIR— > oo. As it can be expected, the proposed and the worst-case-based methods perform almost 
equivalently. 

C. Example 3: Signal Spatial Signature Mismatch Due to Coherent Local Scattering 

The third example corresponds to the scenario of coherent local scattering ||34~|| . In this case, the 
desired signal steering vector is distorted by local scattering effects so that the presumed signal 
steering vector is a plane wave, whereas the actual steering vector is formed by five signal paths 
as 

4 

a = p + £V^b(0 i ) (55) 

i=l 

where p corresponds to the direct path and b(0j), i = 1,2,3,4 correspond to the coherently 
scattered paths. We model the ith path b(0j) as a a plane wave impinging on the array from the 
direction 6^. The angles 9i, i = 1,2,3,4 are independently drawn in each simulation run from a 
uniform random generator with mean 3° and standard deviation 1°. The parameters ipi, i = 1, 2, 3, 4 
represent path phases that are independently and uniformly drawn from the interval [0, 2n] in each 
simulation run. Note that 6i and ipi, i — 1, 2, 3, 4 change from run to run but do not change from 
snapshot to snapshot. 

Fig. [6] displays the performance of all four methods tested versus the number of training snapshots 
K for fixed single-sensor SNR = 20 dB. Note that the SNR in this example is defined by taking 
into account all signal paths. The performance of the same methods versus SNR for the fixed 
training data size K = 30 is displayed in Fig. [7J Similar to the previous example, the proposed 
beamformer significantly outperforms other beamformers due to its ability to estimate the actual 
steering vector with a hight accuracy. 

VII. Conclusion 

A new approach to robust adaptive beamforming in the presence of signal steering vector errors 
has been developed. According to this approach, the actual steering vector is first estimated using 
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its presumed (prior) value, and then this estimate is used to find the optimal beamformer weight 
vector. The problem of signal steering vector estimation belongs to the class of homogeneous 
QCQP problems. It has been shown that this problem can be solved using the SDP relaxation 
technique or the strong duality theory and the exact solution for the signal steering vector can 
be found efficiently or even in some cases in closed-form. As compared to another well-known 
robust adaptive beamforming method based on the signal steering vector correction/estimation, 
that is, the eigespace-based method, the proposed technique does not suffer from the subspace 
swap phenomenon since it does not use eigenvalue decomposition of the sample covariance matrix. 
Moreover, it does not require any knowledge on the number of interferences. As compared to 
the well-known worst-case-based and probabilistically-constrained robust adaptive beamformers, 
the proposed technique does not use any design parameters of non-unique choice. Our simulation 
results demonstrate the superior performance for the proposed method over the aforementioned 
robust adaptive beamforming methods for several frequently encountered types of signal steering 
vector errors. 



Proof of Theorem 1 

Let a be a feasible point for the original problem (fT5l)-(fT7l). It is straightforward to see that 
A = kk H is also a feasible point for the relaxed problem (f25l)— <f28T) . Thus, the necessity statement 
of the theorem follows trivially. 

Let us now prove sufficiency. Let A = YldLi Aibjbf be a feasible point for the relaxed problem 
(f25l)— (f28T) . where Aj and bj, i = 1, • • • ,M are, respectively, the eigenvalues and eigenvectors of A. 
Also let / = argmhij bfCbj and a = \fMhi. Then the following holds 



Appendix 



a"a= Mbf b, 



M 



(56) 



and the constraint (TTBl) is satisfied 



Moreover, we can write that 




(57) 



where J2i=i ^« = Tr(A) = M. Using the following inequality 




(58) 
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and (1571) , we obtain that 

M 

k H Ck< ^A.bfCb,. (59) 
The right hand side of (|59l ) can be further rewritten as 

M M 

Aibf Cbt = \Tr(Cbib?). (60) 

i=l i=l 

Using the property of the trace that a sum of traces is equal to the trace of a sum, we obtain that 

M I K \ 

]T XiTriCbibf ) = Tr C ^(A^bf ) . (61) 

i=l \ i=l / 

Moreover, since 2~2f=i ^jbjbf^ = A, we have 

M / A' \ 

^ Aibf Cb, = Tr Aibibf = Tr(CA). (62) 

i=l \ i=l / 

Substituting (|62l ) in the left hand side of (|59l , we finally obtain that 

a^Ca < Tr(CA) < A . (63) 

Therefore, the constraint (fTTT) is also satisfied and, thus, a = \fMhi is a feasible point for (fT5l)-(fT7l) 
that completes the proof. □ 

Proof of Theorem 2 

Let A* be the optimal minimizer of the relaxed problem (f25l) — (f28l). and its rank be r. Consider 
the following decomposition of A* 

A* = YY H (64) 

where Y is an iV x r complex valued matrix. It is trivial that if the rank of the optimal minimizer 
A* of the relaxed problem (f25l) — (f28l) equals one, then Y is the optimal minimizer of the original 
problem (fT5l)-(fT7l). Thus, it is assumed in the following that r > 1. 
We start by considering the following auxiliary optimization problem 

min Tr(R _1 YGY H ) (65) 

G 

subject to Tr(YGY H ) = M (66) 

Tr(CYGY^) = Tr(CA*) (67) 

G y (68) 
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where G is an r x r Hermitian matrix. The matrix A in (|25l)-(|28T) can be expressed as a function 
of the matrix G in (I65l)-(l68l) as A(G) = YGY H . Moreover, it can be easily shown that if G is a 
positive semi-definite matrix, then A(G) is also a positive semi-definite matrix. In addition, if G 
is a feasible solution of (|65T)-(|68T), A(G) is also a feasible solution of (|25l)-(|28l). The latter is true 
because A(G) is a positive semi-definite matrix and it satisfies the constraints Tr(A(G)) = M and 
7V(C A(G)) = Tr(CA*) < A . This implies that the minimum value of the problem (l65T)-(l68l) 
is greater than or equal to the minimum value of the problem (l25l)-([28T). 

It is then easy to verify that G* = I r is a feasible point of the auxiliary optimization problem 
<[65]>-<|68]>- Moreover, the fact that Tr(R" 1 YG*Y H ) = Tr(R _1 A*) = (3 (here (3 denotes the 
minimum value of the relaxed problem (f25T)— <f28T)^) together with the fact that the minimum value of 
the auxiliary problem (I65l)-(l68l) is greater than or equal to (3, implies that G* = I r is the optimal 
solution of the auxiliary problem (I65l)-(l68l). 

Next, we show that if G' is a feasible solution of (I65l)-(l68l). then it is also an optimum minimizer 
of (I65l)-(l68l). Therefore, A(G') = YG'Y H is also an optimum minimizer of (l25l)-(l28l). Towards 
this end, let us consider the following dual to (I65l)-(l68l) problem 

max u x M + u 2 Tr(CA*) (69) 
subject to Y H R _1 Y - u x Y H Y - u 2 Y H CY y Z (70) 
Z y (71) 

where v x and v 2 are the Lagrange multipliers associated with the constraints (l66l) and (1671) . respec- 
tively, and Z is an r x r Hermitian matrix. Note that the optimization problem (I65l)-(l68l) is convex. 
Moreover, it satisfies the Slater's conditions because, as it was mentioned, the positive definite 
matrix G = I r is a feasible point for (I65l)-(l68l). Thus, the strong duality between (|65T)-(|68T) and 
©-(EB holds. 

Let i/*, u 2 , and Z* be one possible optimal solution of the dual problem (l69l)-(r7TI). Since strong 
duality holds, u^M+u^Tr^A*) = and I r is an optimal solution of the primal problem (I65l)-(l68l). 
Moreover, the complementary slackness condition implies that 

Tr(G*Z*) = Tr(Z*) = 0. (72) 

Since Z* has to be a positive semi-definite matrix, the condition (1721) implies that Z* = 0. Then it 
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follows from (TO} that 

Y H R _1 Y h v* x Y H Y + v* 2 Y H CY. (73) 

Using the fact that Y H Y and Y H CY are positive semi-definite matrices, it can be easily verified 
that the constraint (1731) is active, i.e., it is satisfied as equality for optimal v\ and v\. Therefore, 
we can write that 

Y^R *Y = v\Y H Y + i/*Y H CY. (74) 

Let G' be another feasible solution of (I65l)-(l68l) different from I r . Then the following conditions 
must hold 

Tr(Y H YG') = M (75) 

Tr(Y H CYG') = Tr(CA*) (76) 

G' y 0. (77) 

Multiplying both sides of the equation ( 1741) by G', we obtain 

Y H R- 1 YG' = ulY H YG' + ^Y fl CYG'. (78) 

Moreover, taking the trace of the right hand and left hand sides of (1781) . we have 

rr(Y H R _1 YG') = i/JTr(Y H YG') + ^Tr(Y^CYG') 

= v\M + v* 2 Tr(CA*) = /3. (79) 

This implies that G' is also an optimal solution of the problem (I65l)-(l68l). Therefore, every feasible 
solution of (I65l)-(l68l) is also an optimal solution. 

Finally, we show that there exists a feasible solution of (l65l)-(l68l) whose rank is one. As it 
has been proved above, such feasible solution will also be optimal. Let G = vv^. Thus, we are 
interested in finding such v that 

Tr(Y H Yw H ) = M (80) 
Tr(Y H CYvv H ) = Tr(CA*). (81) 
Equivalently, the conditions (1801) and (ISTI) can be rewritten as 

v^Y^Yv = M (82) 
v^Y^CYv = Tr(CA*). (83) 
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We can further write that 

-v ff Y fl Yv = 1 (84) 
M V ; 

v g Y v = 1. (85) 

Tr(CA*) 

Moreover, equating the left hand side of (f8~4l) to the left hand side of (l85l) . we obtain that 

— v ff Y H Yv = v g Y ^ Y v. (86) 
M Tr(CA*) 

Subtracting the left hand side of < [8"6l ) from its right hand side, we also obtain that 

v^CY \ 

= v = v ff Dv = 0. (87) 

Tr (CA*) J 

Considering the fact that Tr(D) = 0, the vector v can be chosen as the sum of the eigenvectors of 
the matrix D in d29l ). Note that v can be chosen proportional to the sum of the eigenvectors of the 
matrix D such that v^Y^Yv = M is satisfied. It will also imply that v^Y^CYv = Tr(CA*) 
and, thus, (|82l) and (l83~l) are satisfied. 

So far we have found a rank one solution for the auxiliary optimization problem (|66l)-(l68l), 
that is, G = vv^. Since G = vv ff is the optimal solution of the auxiliary problem (I66l)-(l68l). 
then YGY H = (Yv)(Yv) ff is the optimal solution of the relaxed problem (l25l)-(l28l). Moreover, 
since the solution (Yv)(Yv) /f is rank-one, Yv is the optimal solution of the original optimization 
problem (TT3T) — (fTvT) . This completes the proof. □ 

Proof of Theorem 3 

Let A* be one optimal solution of the problem (l25l)-(l28l) whose rank r is greater than one. Using 
the rank-one decomposition of Hermitian matrices 11351 . the matrix A* can be written as 



where 



A* = J2 z j z f ( 88 ) 

3=1 



1„ , . „ M 



zfz j = -Tr(A*) = — , j = 1, . . . , r (89) 

zf C Zj = ^Tr(CA*), j = l,...,r. (90) 

Let us show that the terms z^R~ 1 z J , j = 1, . . . , r are equal to each other for all j = 1, . . . , r. We 
prove it by contradiction assuming first that there exist such z m and z„, m ^ n that z^R _1 z m < 
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z„ R _1 z n . Let the matrix Aq be constructed as Aq = A* — z n z„ + z m z^. It is easy to see that 
Tr(A*) = Tt(Aq) and Tr(CA*) = T^CAq), which means that Aq is also a feasible solution 
of the problem (|25l)-(|28T). However, based on our assumption that z^R -1 z m < z^R _1 z n , it can 
be concluded that Tr (R _1 Aq) < Tt^Rt^Aq) that is obviously a contradiction. Thus, all terms 
z^R~ 1 z J , j — 1, . . . , r must take the same value. Using this fact together with the equations ( |89l ) 
and (1901) , we can conclude that rzjz^ for any j = 1, . . . , r is the optimal solution of the relaxed 
optimization problem (T25l)-(l28l) which has rank one. Thus, the optimal solution of the original 
problem (fT5l)-(fT7l) is -\A" z j f° r an Y j = 1> ■■■,r. Since the vectors Zj, j = l,...,r in (1881) are 
linearly independent and each of them gives an optimal solution to the problem (IT3T) — (ITvT). we 
conclude that the optimal solution to (fT5l)-(fT7l) is not unique up to a phase rotation. However, it 
contradicts the assumption that the optimal solution of (fT3T) — (TTvT) is unique up to a phase rotation. 
Thus, the optimal solution A* to the relaxed problem (l25T)-(l28l) must be rank-one. This completes 
the proof. □ 
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Fig. 1. The term d H (8)Cd(8) in the constraint d 14b versus different angles. 
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Fig. 2. Output SINR versus training sample size K for fixed SNR = 20 dB and INR = 30 dB . 
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Fig. 3. Output SINR versus SNR for training data size of K = 30 and INR = 30 dB. 
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Fig. 4. Output SINR versus training sample size K for fixed SNR = 20 dB and INR = 30 dB. 
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Fig. 5. Output SINR versus SNR for training data size of K = 30 and INR = 30 dB. 



— OPTIMAL SINR 

- - WORST-CASE-BASED BEAMFORMER 

— EIGENSPACE-BASED BEAMFORMER 

- - PROPOSED BEAMFORMER 

SQP-BASED BEAMFORMER 



NUMBER OF SNAPSHOTS 



Fig. 6. Output SINR versus training sample size K for fixed SNR = 20 dB and INR = 30 dB. 
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Fig. 7. Output SINR versus SNR for training data size of K = 30 and INR = 30 dB. 



August 6, 2010 



DRAFT 



